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JUMPING AROUND THE CENTER OF K-SPACE IN ANNULAR 

SEQUENTIAL STEPS 

Background of the Invention 

5 The following relates to the magnetic resonance arts. It finds particular 

application in magnetic resonance imaging during a transient magnetic contrast agent 
bolus, and will be described with particular reference thereto. However, it also finds 
application in magnetic resonance imaging of transient events generally, such as in 
gated magnetic resonance imaging. 
10 In certain magnetic resonance imaging applications, k-space data is preferably 

acquired during a brief time interval. For example, in contrast-enhanced imaging the 
magnetic resonance imaging data are preferably acquired in the typically brief time 
interval during which the magnetic contrast agent bolus provides the most strongly 
enhanced magnetic resonance contrast in the tissue* In some contrast enhanced 
15 magnetic resonance angiography techniques, center of k-space data are collected while 
the contrast agent is in the arteries and before it arrives in the veins to provide good 
visualization of arterial flow by suppressing the not yet enhanced veneous signal. As 
another example, in cardiac gated imaging the magnetic resonance imaging data are 
preferably acquired in a short time interval corresponding to a selected cardiac phase so 
20 as to reduce motion artifacts. 

It will be appreciated that in these and similar applications the time interval for 
acquiring data is brief, but generally not well-defined. A contrast agent bolus, for 
example, typically produces a ramp-up in magnetic resonance signal intensity as the 
contrast agent washes into the imaged region. The slope of the signal intensity rise 
25 begins to level off until eventually a peak magnetic resonance signal intensity is 
reached. Beyond this intensity peak, the signal intensity begins to decrease as the 
contrast agent washes out of the imaged region. Hence, the brief time interval for 
optimal imaging is not well-defined, but it is usually advantageous to concentrate 
acquisition of data around the time of the maximum bolus signal intensity. For 
30 example, separation of arterial and veneous signal is, achieved by acquiring the center 
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of k-space data while the bolus induced signal maximum occurs in the vessel of 
interest Similarly, in gated cardiac imaging, the objective is to image a single cardiac 
phase; as the imaging time increases, the cardiac muscle increasingly deviates from the 
target cardiac phase. 

5 To improve imaging of a transient event, it is known to order the k-space 

sampling so that k-space samples around the cento: of k-space are acquired closest to 
the optimal time, with k-space samples away from k-space center being acquired at less 
optimal times. For example, Riederer et al., U.S. patent no. 5,122,747 discloses a type 
of centric sampling in which the sampling trajectory spirals outward from k-space 

10 center, so that the initial k-space samples acquired at the optimum time are near 
k-space center. The k-space samples close to the center of k-space are low frequency 
samples that contribute substantially to the overall image. By distinction, k-space 
samples at larger distances from the center of k-space are high frequency samples that 
contribute to finer details of the image. Hence, by acquiring samples near the center of 

15 k-space first, the samples that most contribute to the overall image are acquired closer 
to the bolus peak, target arterial or veneous phase, target cardiac phase, or other 

optimal time interval. 

Although data acquisition techniques exist which advantageously concentrate 
data acquisition toward k-space center, these existing techniques have certain 

20 disadvantages. The data are not collected using a geometry that comports with typical 
Fourier transform-based reconstruction algorithms. Hence, a substantial amount of 
sorting of k-space samples is usually required to prepare the collected imaging data for 
image reconstruction processing. Such sorting can be computationally intensive, 
slowing down the data processing and delaying generation of useful images. 

25 Moreover, by collecting data spiraling out from the center, substantially the 

entire data set must be collected before there is a completely acquired plane of k-space 
on which to begin Fourier-based image reconstruction. This prevents concurrent data 
acquisition and image reconstruction processing, which further delays the generation of 
useful images. 

30 Yet another problem with existing techniques is that they typically involve 

applying a substantial number of large magnetic field gradient steps in order to follow 
the complex centric or other selected k-space acquisition trajectory. These large 
magnetic field gradients can generate eddy currents that contribute to image artifacts. 
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The present invention contemplates improved apparatuses and methods that 
overcome the aforementioned limitations and others. 

* 

Brief Summary of the Invention 

According to one aspect, a magnetic resonance imaging method is provided* A 
5 k-space is divided into a central region disposed at k-space center and one or more 
annular surrounding regions having increasing distances from k-space center. The one 
or more annular surrounding regions including an outermost surrounding region having 
a largest distance from k-space center. The k-space samples in the central region are 
acquired. Subsequent to the acquiring of k-space samples in the central region, k-space 
10 samples in the one or more annular surrounding regions (104, 106) are acquired. The 
k-space samples in the outermost surrounding region (106) are acquired last The 
k-space samples are acquired in at least the outermost surrounding region using a 
row-by-row data acquisition ordering in which each row of k-space samples acquired in 
the outermost surrounding region, together with selected already-acquired k-space data 
15 from the regions other than the outermost surrounding region, forms a completed data 
set for reconstructing an image plane. Each completed data set is reconstructed into a 
reconstructed image plane without waiting for all k-space samples in the outermost 
surrounding region to be acquired such that the reconstructing occurs at least partially 
concurrently with the acquiring. 
20 According to another aspect, a magnetic resonance imaging apparatus is 

disclosed. A magnetic resonance imaging scanner images an associated imaging 
subject A magnetic resonance imaging controller performing a method including: (i) 
dividing k-space into a central region disposed at k-space center and one or more 
annular surrounding regions having increasing distances from k-space center and 
25 including an outermost surrounding region of largest distance from k-space center; (ii) 
determining an optimum time for imaging a magnetic contrast agent bolus; (iii) 
acquiring k-space samples in the central region at about the optimum time, and (iv) 
after acquiring the k-space samples in the central region, acquiring k-space samples in 
the one or more annular surrounding regions. The acquiring in at least the outermost 
30 surrounding region uses a plane-by-plane data acquisition ordering in which all k-space 
samples in the outermost surrounding region belonging to a current k-space plane are 
acquired to complete the current k-space plane before samples in the outermost 
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surrounding region belonging to other k-space planes are acquired. A reconstruction 
processor reconstructs the completed current k-space plane into a reconstructed plane 
image without waiting for other k-space planes to be completed. 

According to yet another aspect, a magnetic resonance imaging apparatus is 

5 disclosed. A means is provided for dividing k-space into a central region disposed at 
k-space center and one or more annular surrounding regions having increasing 
distances from k-space center. The one or more annular surrounding regions include an 
outermost surrounding region having a largest distance from k-space center. A means is 
provided for acquiring k-space samples in the k-space. The k-space samples in the 

10 central region are acquired first. The k-space samples in the outermost surrounding 
region are acquired last The k-space samples in at least the outermost surrounding 
region are acquired using a row-by-row data acquisition ordering in which each row of 
k-space samples acquired in the outermost surrounding region completes a k-space 
plane. A means is provided for reconstructing each completed k-space plane into a 

15 reconstructed image plane without waiting for all k-space samples in the outermost 
surrounding region to be acquired. 

One advantage resides in shortened imaging session times. 
Another advantage resides in reduced image artifacts. 

Another advantage resides in enabling concurrent imaging data acquisition, 
20 reconstruction, and display processing. 

Yet another advantage resides in reduced data sorting complexity. 

Still yet another advantage resides in reduced generation of eddy currents. 

Numerous additional advantages and benefits will become apparent to those of 
ordinary skill in the art upon reading the following detailed description of the preferred 
25 embodiments. 

Brief Description of the Drawings 
The invention may take form in various components and arrangements of 
components, and in various process operations and arrangements of process operations. 
The drawings are only for the purpose of illustrating preferred embodiments and are 
30 not to be construed as limiting the invention. 
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FIGURE 1 diagrammatically shows a magnetic resonance imaging system. 

FIGURE 2 shows the division of the k-space sampling into (i) a central region 
disposed at k-space center, and (ii) two annular surrounding regions. 

FIGURE 3 shows a suitable random or pseudorandom acquisition ordering for 
5 k-space samples lying in the central k-space region. 

FIGURE 4 shows a suitable row-by-row acquisition ordering for k-space 
samples lying in the inner annular surrounding k-space region. 

FIGURE 5 shows a suitable row-by-row acquisition ordering for k-space 
samples lying in the outermost annular k-space region. 
10 FIGURE 6 diagrammatically shows relative timing of the data acquisition, the 

data reconstruction, and the data display aspects of the imaging. 

FIGURE 7 shows a suitable embodiment of the reconstruction processor of 
FIGURE 1. 

FIGURE 8 shows a suitable random or pseudorandom acquisition ordering for 
1 5 k-space samples lying in the inner annular surrounding k-space region. 

, Detailed Description of the Preferred P mhftrifrnents 
With reference to FIGURE 1, a magnetic resonance imaging scanner 10 
includes a housing 12 defining a generally cylindrical scanner bore 14 inside of which 
an associated imaging subject 16 is disposed. Main magnetic field coils 20 are disposed 
inside the housing 12, and produce a main Bo magnetic field directed generally parallel 

■ 

20 with a central axis 22 of the scanner bore 14. The main magnetic field coils 20 are 
typically superconducting coils disposed inside cryoshrouding 24, although resistive 
main magnets can also be used. The housing 12 also houses or supports magnetic field 
gradient coils 30 for selectively producing magnetic field gradients in the bore 14. The 
housing 12 further houses or supports a radio frequency body coil 32 for selectively 

25 exciting and/or detecting magnetic resonances. The housing 12 typically includes a 
cosmetic inner liner 36 defining the scanner bore 14. For some imaging applications, 
one or more local radio frequency coils such as a surface coil 44 may be used for radio 
frequency excitation and/or detection. 

In contrast-enhanced magnetic resonance imaging, a bolus injector 46 is used to 

30 rapidly infuse a selected magnetic contrast agent, such as a gadolinium-based contrast 
agent, into the imaging subject 16. For a bolus injection into the bloodstream of a 
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human imaging subject, the contrast agent rapidly washes into and out of tissue of 
interest. The bolus transit time between wash-in and wash-out is typically on the order 
of a few seconds. In the illustrated embodiment, the bolus injector 46 is coupled into 
the bloodstream of the imaging subject 16 by a catheter 48. 
5 A magnetic resonance imaging controller 50 operates magnetic field gradient 

controllers 52 to selectively energize the magnetic field gradient coils 30, and operates 
a radio frequency transmitter 54 coupled to the radio frequency coil 32 as shown, or 
coupled to a local coil such as the surface coil 44, to selectively inject radio frequency 
excitation pulses into the imaging subject 16. By selectively operating the magnetic 
10 field gradient coils 30 and the radio frequency coil 32, magnetic resonance is generated 
and spatially encoded in at least a portion of a region of interest of the imaging subject 
16. During imaging data acquisition, the magnetic resonance imaging controller 50 
operates a radio frequency receiver 56 coupled to the radio frequency surface coil 44 as 
shown, or coupled to the whole body coil 32, to acquire magnetic resonance samples 
1 5 that are stored in a magnetic resonance data buffer 60. 

The acquired imaging data are reconstructed by a reconstruction processor 62 
into an image representation, typically using a Fourier transform-based reconstruction 
algorithm. The reconstructed image or images generated by the reconstruction 
processor 62 are stored in an images memory 64, and can be displayed on a display 
20 device 66 of a user interface 68, stored in non-volatile memory, transmitted over a local 
intranet or the Internet, viewed, stored, manipulated, or so forth. The user interface 68 
also preferably enables a radiologist, technician, or other operator of the magnetic 
resonance imaging scanner 10 to communicate with the magnetic resonance imaging 
controller 50 to select, modify, and execute magnetic resonance imaging sequences. 

* 

25 The described magnetic resonance imaging system is an illustrative example 

only. The magnetic resonance data acquisition techniques described herein can be used 
with substantially any type of magnetic resonance imaging scanner, including but not 
limited to horizontal bore scanners, vertical bore scanners, open magnet scanners, and 
so forth. Moreover, the magnetic resonance data acquisition techniques described 

30 herein can be employed in conjunction with arrays of receive coils for parallel imaging, 
SENSE imaging, and the like. 

In contrast enhanced imaging, the bolus injector 46 administers a magnetic 
contrast agent bolus to the imaging subject 16. When administered into the 
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bloodstream, the bolus washes into the imaging region of interest, and then washes out, 
typically on the order of a few seconds. Accordingly, the imaging data acquisition is 
preferably timed to acquire imaging data at the peak of the magnetic contrast provided 
by the bolus in the imaging region of interest As will be described, the imaging data is 
5 acquired in a manner in which k-space samples in a region near the center of k-space 
are acquired first, followed by samples in regions at larger distances from k-space 
center. In one approach for timing the imaging, the magnetic resonance imaging 
controller 50 sends a bolus triggering signal 70 to the bolus injector 46. The bolus 
injector 46 administers the bolus via the catheter 48 responsive to receiving the bolus 
10 triggering signal 70. The magnetic resonance imaging controller 50 initiates die 
magnetic resonance imaging data acquisition at a pre-selected time after sending the 
bolus triggering signal 70, where the pre-selected time corresponds to die time it takes 
for the administered bolus to wash into the imaging region of interest 

In another approach for timing the imaging data acquisition, the scanner 10 
15 excites and acquires a projection or other magnetic resonance signal from the region of 
interest The magnetic resonance imaging controller 50 monitors the intensity of this 
signal. The bolus injector 46 administers the magnetic contrast agent bolus. As the 
administered bolus washes into the imaging region of interest, the intensity of the 
projection or other magnetic resonance signal changes accordingly. The magnetic 
20 resonance imaging controller 50 detects this change and begins image data acquisition. 
In yet other embodiments, such as in cardiac gated imaging, the magnetic resonance 
imaging controller 50 uses a suitable cardiac gating signal, for instance extracted from 
an electrocardiograph signal, to trigger data acquisition. Similarly, substantially any 
other physiological event can be detected and used to trigger initiation of magnetic 
25 resonance imaging of a volume of interest affected by the physiological event 

With continuing reference to FIGURE 1 and with further to FIGURES 2-5, the 
magnetic resonance imaging controller 50 divides k-space 100 into a central region 102 
disposed at k-space center and one or more annular surrounding regions having 
increasing distances from k-space center, namely two annular surrounding regions 104, 
30 106 in the example embodiment illustrated in FIGURE 2. The k-space representations 
of FIGURES 2-5 show a view perpendicular to a primary data acquisition coordinate 
designated n Cp n , and perpendicular to an orthogonal secondary data acquisition 
coordinate designated "Cs". An unlabeled readout coordinate is oriented transverse to 
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illustraled the primary Cp and secondary Cs coordinates, and hence parallel with the 
direction of view. In some embodiments, the Cp coordinate is a slice-select coordinate 
while the Cs coordinate is a phase-encode coordinate; however, the primary and 
secondary coordinates can have substantially arbitrary spatial orientations. For 
5 example, the phase-encode direction could be designated the primary coordinate Cp 
while the slice-select direction could be designated the secondary coordinate Cs, or one 
or both of primary Cp and secondary Cs coordinates may be inclined to the conventional 
sagittal, coronal, and axial anatomical coordinates. Moreover, while in the illustrated 
embodiment the primary Cp, secondary c s , and readout data acquisition coordinates are 
10 mutually orthogonal, in some contemplated embodiments a non-orthogonal data 
acquisition coordinates system may be employed. 

In FIGURES 2-5, a rectangular grid of k-space samples 110 is illustrated. In 
these FIGURES, each k-space sample or readout line is depicted as a dot. The 
illustrated grid is 10*10 for simplicity of illustration. In practice, grids of 128*128 to 
15 1024*1024 are routine. In embodiments in which a three-dimensional volume is 
imaged, these k-space samples 110 each correspond to a readout line of k-space data 
mat extends along the readout coordinate (mat is, extends transverse to the primary Cp 
and secondary Cg coordinates, and thus appears as a dot from the perspective of 
FIGURES 2-5). Although a square 10*10 grid of samples is illustrated, the number of 
20 samples in the Cp and Cs directions can be different In such three-dimensional volume 
imaging, the central region 102 will extend into the readout direction to define a 
generally cylindrical central region, while the annular regions 104, 106 will similarly 
extend into the readout direction to define three-dimensional annular volumes 
surrounding the cylindrical central region volume 102. 
25 The magnetic resonance imaging controller 50 acquires the central region 102 

first using an acquisition ordering shown in FIGURE 3, and then acquires the inner 
annular surrounding region 104 using an acquisition ordering shown in FIGURE 4, and 
finally acquires die outermost annular surrounding region 106 using an acquisition 
ordering shown in FIGURE 5. Some preferred sampling orderings for use in acquiring 
30 data from each of these regions 102, 104, 106 is next described. 

With particular reference to FIGURE 3, each k-space sample 110 in the central 
region 102 is labeled with its ordinal number in the acquisition sequence. The example 
central region 102 encompasses twelve k-space samples labeled with ordinal numbers 
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1, 2, 3, ... 12. In the illustrated embodiment, a substantially random or pseudorandom 
ordering of the k-space samples acquisition is used for acquiring the k-space samples in 
die central region 102. A random or pseudorandom ordering for acquiring samples of 
the central region 102 advantageously provides robustness with respect to errors in 
5 triggering the data acquisition. If the magnetic resonance data acquisition is triggered 
slightly early or slightly late such that the peak intensity contrast provided by the 

* 

magnetic contrast agent bolus does not precisely coincide with acquisition of the 
central region 102, then this random or pseudorandom ordering reduces the sensitivity 
of the image quality to this triggering error. In this context, "random or pseudorandom" 

10 is intended to encompass, in addition to the ordinary meaning of those terms, any 
pre-determined ordering selected to have statistical characteristics similar to those of a 
random or pseudorandom ordering. 

A random or pseudorandom ordering typically has some cost in terms of 
increased magnetic field gradients employed to follow large jumps between the random 

15 sampling positions. However, because the central region 102 is generally relatively 
small and close to the center of k-space, the magnetic field gradients required to make 
die random jumps within the central region 102 are typically acceptable. Alternatively, 
a centric or other ordering that collects the k-space samples closest to die center of 
k-space can be used in die central region 102 instead of the illustrated random or 

« 

20 pseudorandom ordering. 

■ 

With particular reference to FIGURE 4, once all the k-space samples in the 
central region 102 are acquired, the magnetic resonance imaging controller 50 acquires 
k-space samples in the inner surrounding annular region 104 using a serpentine 
row-by-row acquisition ordering that is designed to acquire a row of readout samples 

25 having a given Cp coordinate value as a group (that is, designed to acquire a row of 
readout samples along die Cs coordinate direction) while limiting die magnitudes of the 
applied magnetic field gradients. The acquisition order of the k-space samples is 
represented in FIGURE 4 where each acquired row of samples is indicated by an arrow 
labeled by its ordinal number of plane acquisition. In FIGURE 4, there are six rows 

30 labeled 1, 2, 3, 4, 5, 6. The direction of the arrow indicates the direction the row of 
samples is traversed during data acquisition. The first acquired sample is at the bottom 
of FIGURE 4 around the center position of the primary Cp coordinate. Starting at this 
k-space sample, acquisition of a row of samples traverses a positive or upward 
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direction in c* coordinate, and does not include acquisition of those k-space samples 
along the row which lie inside the central region 102 and which thus have already been 
required It will be appreciated that traversing this row of k-space samples involves 
gradient steps in the secondary Cs coordinate direction, without application of any 

5 gradients in the primary Cp direction. 

Once the first row of k-space samples are acquired, a gradient pulse in the 
primary Cp direction is applied to move to another row, which also lies near the center 
position of the primary Cp coordinate. This second row of k-space samples is traversed 
in the negative or downward direction, that is, the second row is traversed in the 

10 opposite direction from the first row. The second row is also traversed using only 
secondary ^ coordinate gradient pulses, without application of any gradients in the 

primary Cp direction. 

Once the second row k-space samples is acquired, a gradient pulse in the 
primary Cp direction is applied to move to the third row, which lies further out from the 
15 center position of the primary Cp coordinate. The third row of samples is traversed in 
the positive or upward direction. Once the third row of k-space samples is acquired, a 
gradient pulse in the primary Cp direction is applied to move to the fourth row, which is 
traversed in the negative or downward direction. Once the fourth row of k-space 
samples is acquired, a gradient pulse in the primary Cp direction is applied to move to 
20 die fifth row, which lies at an edge of the surrounding annular region 104. The fifth row 
of k-space samples is traversed in the positive or upward direction. Once the fifth row 
of k-space samples is acquired, a gradient pulse in the primary Cp direction is applied to 
move to the sixth and final row of k-space samples intersecting the surrounding region 
104, which lies at the opposite edge of the surrounding annular region 104. The sixth 
25 row of k-space samples is traversed in the negative or downward direction. 

The described row-by-row acquisition pattern in the inner surrounding annular 
region 104 has certain advantages. In three-dimensional Fourier-encoded imaging, each 
k-space sample 110 corresponds to a readout line of k-space in the readout direction. 
Hence, a row of k-space samples corresponds to a plane of k-space. This plane of 
30 k-space can be transformed into a reconstructed image plane using two-dimensional 
Fourier-based reconstruction processing, thus eliminating data processing delays 
caused by complex sorting operations. Moreover, the row-by-row acquisition ordering 
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limits the magnitudes of changes in amplitude of the magnetic field gradients used in 
traversing the trajectory, which reduces eddy currents that can lead to image artifacts. 

With particular reference to FIGURE 5, once all the k-space samples in the 
central region 102 and in the inner annular surrounding region 104 are acquired, the 
5 magnetic resonance imaging controller SO acquires k-space samples in the outermost 
surrounding annular region 106 using the same type of row-by-row acquisition 
ordering used in acquiring samples in the inner region 104. The acquisition of each row 
of k-space samples is represented in FIGURE 5 by an arrow labeled by its ordinal 
number of acquisition. There are ten row of k-space samples labeled 1, 2, 3, . .., 10. The 
10 direction of the arrow indicates the direction the row of k-space samples is traversed 
during the data acquisition. In similar fashion to the acquisition ordering used in the 
inner surrounding region 104, the row of the outermost surrounding region 106 are 
acquired such that the first and second rows having primary Cp coordinate values 
closest to the center of the outermost region 106 are acquired first, and rows with 
15 increasing primary coordinate Cp value are subsequently acquired, finishing with 
acquisition of the ninth and tenth rows of k-space samples located at the opposite edges 
of the outermost annular region 106, respectively. The direction of acquisition of 
samples along a row in the secondary direction reverses at the beginning of each new 
row to limit applied magnetic field gradient magnitudes. Tims, rows 1, 3, 5, 7, 9 are 
20 traversed in the positive or upward direction, while rows 2, 4, 6, 8, 10 are traversed in 
the negative or downward direction. When acquiring the k-space data of the outermost 
surrounding region 106, k-space samples along the row which lie in the central region 
102 or in the inner annular surrounding region 104, and which have thus already been 
acquired, are not re-acquired during the row-by-row sampling of the outermost 
25 surrounding region 106. 

With continuing reference to FIGURES 1-5 and with further reference to 
FIGURE 6, the row-by-row acquisition of the outermost annular surrounding region 
106 advantageously enables image reconstruction to initiate before all k-space data is 
acquired, and allows image reconstruction to be performed substantially concurrently 
30 with data acquisition once a full row of k-space samples (corresponding to a plane of 
k-space) has been acquired. As diagrammatically shown in FIGURE 6, the k-space 
samples in the central region 102 are acquired during a time interval 120, followed by 
acquisition of the k-space samples in the inner annular surrounding region 104 during a 
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subsequent time interval 122. At the end of the second acquisition time interval 122, no 
plane of k-space has been fully sampled, since the outermost annular region 106 
includes at least some k-space samples from every plane of k-space in the imaging 
volume. 

5 The row-by-row sampling of the outermost annular region 106 is performed 

during a time interval 124 following the time interval 122. Once the first row of 
k-space samples is acquired, all the data for a full plane of k-space are available, and 
reconstruction of that plane can,be initiated by the reconstruction processor 62. 

By way of example, when the primary Cp direction is the slice-select direction, 
10 the secondary Cs direction is the phase-encode direction, and the unlabeled direction 
perpendicular to Cp and Cs is the readout direction, acquisition of a row of k-space 
samples corresponds to acquiring all phase-encode steps for a selected slice. As soon as 
all of the k-space samples along the first row of the outermost region 106 are acquired 
by traversing the phase-encode steps along the row for that selected slice (jumping over 
15 those phase-encode steps already acquired during acquisition of regions 102, 104), the 
central slice can be reconstructed. When all of the data along row 2 is collected, the 
adjoining slice can be reconstructed, and so forth. Moreover, the same principle applies 
for other two-dimensional planes through three-dimensional k-space. 

Hence, a reconstruction time interval 130 initiates subsequent to completion of 
20 sampling of the first row of k-space samples in the outermost annular region 106. The 
reconstruction may extend beyond the time interval 124 during which the samples in 
the outermost annular region 106 are acquired (as indicated by the curved right-side 
"break" of the reconstruction time interval 130 in FIGURE 6); however, substantial 
data processing speedup is obtained through the partial overlap of the acquisition time 
25 interval 124 and die reconstruction time interval 130. 

With brief returning reference to FIGURE 1, in one suitable embodiment, the 

■ 

magnetic resonance imaging controller 50 sends a reconstruction trigger signal 72 to 
the reconstruction processor 62 each time a row of k-space samples has been 
completely acquired. Receipt of the reconstruction trigger signal 72 causes the 
30 reconstruction processor 62 to begin reconstruction of the corresponding image plane. 
In another approach, the reconstruction processor 62 monitors the sampling data buffer 
60 to determine when acquisition of a row of k-space samples has been completed. 
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With returning reference to FIGURE 6, the data processing does not involve a 
substantial amount of data sorting. The surrounding regions 104, 106 are acquired on a 
row-by-row basis that comports with two-dimensional Fourier reconstruction, so that 
substantial sorting is not required for those data* The data from the central region 102 is 

5 randomized or pseudorandomized and it must be sorted; however, this sorting is 
performed contemporaneous with the acquisition of the data, with each sample placed 
into a position in the data buffer corresponding to its k-space coordinates as the next 
sample is being acquired 

Still further, since the reconstruction processor 62 reconstructs each plane as its 

1 0 k-space data acquisition is completed during the acquisition interval 124, the image of 
that reconstructed plane can be displayed by the display device 66 in a displaying time 
interval 136 that at least partially overlaps the acquisition time interval 124 and the 
reconstruction time interval 130. Thus, the image begins to be displayed even while 
data acquisition is ongoing. 

15 With reference to FIGURE 7, a suitable embodiment of the reconstruction 

processor 62 is described. In some embodiments, the sorting of the data is determined 
prior to the acquisition and the expected acquisition order is provided to the 
reconstruction processor 62 in an array of k-space points determined by the magnetic 
resonance imaging controller 50. Hence, no sorting is performed during reconstruction 

20 beyond the initial placement of each data sample into a location in the data buffer 
corresponding to its position in k-space by means of a data sorter 150. With reference 
to FIGURE 6, this sorting by the reconstruction processor 62 during the time interval 
132 is performed concurrently with the acquiring and reconstructing of the data. 

In three-dimensional imaging, each k-space sample 110 corresponds to a 

25 readout line in the readout direction which is orthogonal to the primary Cp and 
secondary Cs directions. Accordingly, a one-dimensional Fourier transform 152 is 
applied to each k-space sample to recover the spatial content in the readout direction. 

Within the reconstruction processor 62, a slice checker 160 determines when a 
complete row of k-space samples, corresponding to a complete data set for 

30 reconstructing an image plane, has been acquired. Alternately, a reconstruction trigger 
signal 72 (see FIGURE 1) could indicate when an image plane data set has been 
acquired. The completed data set, which has already been sorted by the data sorter 150 
and Fourier transformed in the readout direction by the Fourier transform processor 



-14- 

152, is processed by a second one-dimensional Fourier transform processor 162 
operating in the secondary Cs coordinate direction to recover the spatial content along 
the secondary Cs coordinate direction. The output of the second one-dimensional 
Fourier transform processor 162 is a reconstructed image of the just-completed 

5 primary Cp plane. As acquisition of each row of k-space samples corresponding to an 
image plane is completed during the time interval 124, that just-completed image plane 
data set is similarly reconstructed. 

The dividing up and sampling of k-space illustrated in FIGURES 2-5 is only an 
example. In some contemplated embodiments, the inner surrounding region 104 is 

10 omitted, i.e., the data are divided between a randomly acquired central region and a 
peripheral region that is acquired row-by-row. In other contemplated embodiments, Ihe 
number of surrounding annular regions is greater than two. As shown in FIGURES 2 
and 5, the example outermost annular region 106 has a circular or oval inner boundary 
and a square or rectangular outer boundary (in the Cp-Cg plane). Thus, it will be 

IS appreciated that the annular surrounding regions can have various cross-sectional 
shapes. 

The advantageous at least partial concurrency of data acquisition and 
reconstruction; (and optionally also concurrent displaying) is obtained only during 
acquisition of k-space samples in the outermost region 106, that is, only during the time 

20 interval 124. Accordingly, there are advantages to making the outermost region large so 
that it contains a large fraction of the k-space samples compared with inner regions 
102, 104, thus enabling a large overlap of data acquisition and image reconstruction. 

On the other hand, the row-by-row acquisition employed in the outermost 
region 106 is not strongly biased toward acquiring data near the center of k-space, and 

25 so making the outermost region 106 too large can increase temporal artifacts. To reduce 
temporal artifacts, the central region 102 should be made larger so that central k-space 
samples that contribute strongly to the overall image are acquired first 

As the central region 102 is enlarged, however, the sorting time prior to data 
acquisition is increased, causing an increase in acquisition planning time. Moreover, if 

30 the central region 102 is too large it will contain k-space samples relatively distant from 
k-space center. Including one or more inner surrounding annular regions such as the 
inner surrounding annular region 104, sampled on a row-by-row basis as shown in 
FIGURE 4, provides a balance by biasing the early portion of the acquisition toward 
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k-space center while providing a row-by-row acquisition ordering that is conducive to 

rapid data processing. 

In some embodiments, the size of the central region 102 is determined as 
follows. The magnetic resonance imaging parameters for the acquiring of k-space 
5 samples is determined. These parameters include a data acquisition rate. The time 
interval 120 for acquiring samples in the central region 102 is also selected, for 
example based on an expected duration of the magnetic contrast agent bolus in the 
tissue of interest Based on the data acquisition rate and the selected time interval 120, 
the number of k-space samples of the central region 102 is determined, and a suitable 
1 0 central region 102 corresponding to this number of samples is selected. 

With reference to FIGURE 8, in some embodiments it is contemplated to use an 
ordering other than serpentine row-by-row ordering for the surrounding annular regions 
other than the outermost surrounding annular region. FIGURE 8, for example, shows 
an example of a random or pseudorandom sampling ordering for the inner surrounding 
15 annular region 104. Substituting the random or pseudorandom sampling of FIGURE 8 
for the row-by-row ordering of FIGURE 4 can be advantageous for certain imaging 
tasks. For example, this approach provides additional robustness against triggering time 
errors by extending the range of the random or pseudorandom sampling further out 
away from k-space center, while the separate random or pseudorandom sampling of 
20 each of the central region 102 and the inner surrounding region 104 ensures that the 
random or pseudorandom sampling covers all the k-space samples within the boundary 
of the central region 102 before sampling outside of central region 102. However, it 
will be recognized that the random or pseudorandom sampling of the inner surrounding 
region 104 increases the required data sorting time. 
25 Although not illustrated, in some embodiments it is contemplated to employ a 

row-by-row acquisition ordering for the central region 102 and for the one or more 
surrounding annular regions 104, 106. Row-by-row ordering does not by itself strongly 
bias die early sampling toward k-space center, however, by using a relatively small 
. central region 102 and one or more relatively small inner surrounding annular regions 
30 104 biasing of the early sampling toward k-space center can be achieved. Moreover, 
since in these embodiments all regions including the central region 102 are sampled on 
a row-by-row basis, sorting delays in the subsequent data processing are substantially 
reduced. 
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The invention has been described with reference to the preferred embodiments. 
Obviously, modifications and alterations will occur to others upon reading and 
understanding the preceding detailed description* It is intended that the invention be 
construed as including all such modifications and alterations insofar as they come 
within die scope of the appended claims or die equivalents thereof. 
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Claims 



Having described the preferred embodiments, the invention is now claimed to 

be: 

1. A magnetic resonance imaging method comprising: 

dividing k-space (100) into a central region (102) disposed at k-space center and 
one or more annular surrounding regions (104, 106) having increasing distances from 
k-space center, the one or more annular surrounding regions including an outermost 
surrounding region (106) having a largest distance from k-space center, 

acquiring k-space samples in the central region (102); 

subsequent to the acquiring of k-space samples in the central region, acquiring 
k-space samples in the one or more annular surrounding regions (104, 106), the k-space 
samples in the outermost surrounding region (106) being acquired last, the acquiring of 
k-space samples in at least the outermost surrounding region using a row-by-row data 
acquisition ordering in which each row of k-space samples acquired in the outermost 
surrounding region, together with selected already-acquired k-space data from the 
regions other than the outermost surrounding region, forms a completed data set for 
reconstructing an image plane; and 

reconstructing each completed data set into a reconstructed image plane without 
waiting for all k-space samples in the outermost surrounding region to be acquired such 
that the reconstructing occurs at least partially concurrently with the acquiring. 

2. The method as set forth in claim 1, further comprising: 

displaying each reconstructed image plane once it is available without waiting 
for the reconstructing of other image planes. 

3. The method as set forth in claim 1, further comprising: 

synchronizing the acquiring of k-space samples in the central region (102) with 
a trigger signal, the trigger signal being one of: (i) a selected duration after 
administering a magnetic contrast agent bolus, (ii) detecting a change in a magnetic 
resonance signal intensity due to wash-in of a magnetic contrast agent bolus, (iii) 
detecting a gating signal, and (iv) detecting a selected physiological event. 
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4. Hie method as set forth in claim 1, further comprising: 

selecting a plurality of magnetic resonance imaging parameters for the 
acquiring of k-space samples in the central region (102) and in the one or more annular 
surrounding regions (104, 106), the plurality of magnetic resonance imaging 
parameters including at least a data acquisition rate; and 

determining the central (102) region using (i) the selected plurality of magnetic 
resonance imaging parameters and (ii) a time interval (120) for the acquiring of k-space 
samples in the central region (102). 

5. The method as set forth in claim 1, wherein the central region (102) has a 
round or oval perimeter, and the outermost surrounding region (106) has a round or 
oval inner perimeter and a square or rectangular outer perimeter. 

* 

6. The method as set forth in claim 1, wherein the acquiring of k-space samples 
in the central region (102) uses an acquisition ordering other than a row-by-row 
acquisition ordering. 

7. The method as set forth in claim 6, wherein the acquiring of k-space samples 
in the central region (102) uses a random or pseudorandom ordering. 

8. The method as set forth in claim 7, further comprising: 

synchronizing the acquiring of k-space samples in the central region (102) with 
the administering of a magnetic contrast agent bolus. 

9. The method as set forth in claim 6, further comprising: 

sorting the k-space samples of the central region into a row-by-row ordering. 

10. The method as set forth in claim 1, wherein the one or more annular 
surrounding regions (104, 106) include at least two surrounding regions, and the 
acquiring of k-space samples in the one or more annular surrounding regions other than 
the outermost surrounding region uses a random or pseudorandom ordering. 
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11. The method as set forth in claim 1, wherein the one or more annular 
surrounding regions (104, 106) include at least two surrounding regions, and the 
acquiring of k-space samples in every annular surrounding region including the 
outermost surrounding region uses a row-by-row acquisition ordering. 

« 

12. The method as set forth in claim 1, wherein each k-space sample (100) is a 
readout line of k-space. ■ 

13. The method as set forth in claim 1, wherein the acquiring of k-space 
samples in at least the outermost surrounding region using a row-by-row acquisition 
ordering includes: 

acquiring the k-space samples using a serpentine row-by-row acquisition 
ordering. 

14. The method as set forth in claim 1, wherein the acquiring of k-space 
samples in at least the outermost surrounding region using a row-by-row acquisition 
ordering includes: 

applying secondary coordinate magnetic field gradients to traverse each row of 

k-space samples; and 

switching to each new row of k-space samples by applying a primary coordinate 
magnetic field gradient, the primary coordinate being generally transverse to the 
secondary coordinate. 

15. The method as set forth in claim 1, wherein the acquiring of k-space 
samples in at least the outermost surrounding region using a row-by-row acquisition 
ordering includes: 

(i) acquiring a first row of k-space samples by traversing secondary coordinate 
positions in a positive direction at a first primary coordinate position; 

(ii) applying a primary coordinate magnetic field gradient to move to a second 

primary coordinate position; 

(iii) acquiring the second row of k-space samples by traversing secondary 
coordinate positions in a negative direction at the second primary coordinate position; 
and 
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(iv) repeating (i), (ii), and (iii) to acquire a plurality of rows of k-space samples 
indexed by the primary coordinate. 

16. The method as set forth in claim 15, wherein the primary coordinate is a 
slice coordinate, the secondary coordinate is a phase-encode coordinate orthogonal to 
the slice coordinate, and each k-space sample (100) is a readout line along a third 
coordinate orthogonal to both the slice and phase-encode coordinates. 

17. The method as set forth in claim 1, wherein the acquiring of k-space 
samples in at least the outermost surrounding region using a row-by-row acquisition 
ordering includes: 

(i) acquiring a first contiguous portion k-space samples along a row within the 
outermost annular surrounding region; 

(ii) skipping at least samples along the row contained in the central region; 

(iii) acquiring a second contiguous portion k-space samples along the row 
within the outermost annular surrounding region, the second contiguous portion 
k-space samples along the row being separated from the first contiguous portion 
k-space samples along the row by at least the central region; and 

repeating (i), (ii), and (iii) for each row of the row-by-row acquisition. 

18. A magnetic resonance imaging apparatus comprising: 

a magnetic resonance imaging scanner (10) imaging an associated imaging 
subject (16); 

a magnetic resonance imaging controller (50) performing a method including: 

(i) dividing k-space (100) into a central region (102) disposed at 
k-space center and one or more annular surrounding regions (104, 106) 
having increasing distances from k-space center and including an 
outermost surrounding region (106) of largest distance from k-space 
center, 

(ii) determining an optimum time for imaging a magnetic 
contrast agent bolus, 

(iii) acquiring k-space samples in the central region (102) at 
about the optimum time, and 
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(iv) after acquiring the k-space samples in the central region, 
acquiring k-space samples in the one or more annular surrounding 
regions (104, 106), the acquiring in at least the outermost surrounding 
region (106) using a plane-by-plane data acquisition ordering in which 
all k-space samples in the outermost surrounding region belonging to a 
current k-space plane are acquired to complete the current k-space plane 
before samples in the outermost surrounding region belonging to other 
k-space planes are acquired; and 

a reconstruction processor (62) that reconstructs the completed current 
k-space plane into a reconstructed plane image without waiting for other 
k-space planes to be completed. 

19. The magnetic resonance imaging apparatus as set forth in claim 18, 
wherein the acquiring of k-space samples in the central region (102) includes: 

acquiring k-space samples in the central region (102) using an ordering other 
than a plane-by-plane ordering. 

20. The magnetic resonance imaging apparatus as set forth in claim 18, 
wherein the acquiring of k-space samples in the central region (102) includes: 

acquiring k-space samples in the central region (102) using a random or 
pseudorandom ordering. 

21. The magnetic resonance imaging apparatus as set forth in claim 20, further 
comprising: 

a display device (68) that displays each reconstructed plane image once the 
reconstruction processor (62) completes reconstruction of the corresponding k-space 
plane, without waiting for the reconstruction processor to reconstruct other k-space 
planes. 

22. A magnetic resonance imaging apparatus comprising: 

means for dividing k-space (100) into a central region (102) disposed at k-space 
center and one or more annular surrounding regions (104, 106) having increasing 
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distances from k-space center* the one or more annular surrounding regions including 
an outermost surrounding region (106) having a largest distance from k-space center, 

means for acquiring k-space samples in the k-space (100), the k-space samples 
in the central region (102) being acquired first, the k-space samples in the outermost 
surrounding region (106) being acquired last, the k-space samples in at least the 
outermost surrounding region being acquired using a row-by-row data acquisition 
ordering in which each row of k-space samples acquired in the outermost surrounding 
region completes a k-space plane; and 

means (62) for reconstructing each completed k-space plane into a 
reconstructed image plane without waiting for all k-space samples in the outermost 
surrounding region (106) to be acquired. 

23. The apparatus as set forth in claim 22, wherein the acquiring means 
includes: 

(i) means 1 for acquiring a row of k-space samples by traversing secondary 
coordinate (Cg) positions in a positive direction at a first primary coordinate position 

(cp); 

(ii) means for applying a primary coordinate (cp) magnetic field gradient to 
move to a second primary coordinate (Cp) position; 

(iii) means for acquiring the second row of k-space samples by traversing 
secondary coordinate (c 8 ) positions in a negative direction at the second primary 
coordinate (cp) position; and 

(iv) means for repeatedly invoking the means (i), (ii), and (iii) to acquire a 
plurality of rows of k-space samples indexed by the primary coordinate (cp). 



24. The apparatus as set forth in claim 22, wherein the acquiring means 
acquires at least the central region (102) using other than a row-by-row acquisition 
ordering, and the reconstructing means (62) includes: 

means (150) for sorting k-space samples acquired in other than a row-by-row 
acquisition ordering into a row-by-row ordering. 
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25. The apparatus as set forth in claim 22, wherein the reconstructing means 
(62) includes: 

means (150) for organizing k-space samples of the completed k-space plane 
from the central region (102) and from the one or more annular surrounding regions 
(104, 106) into a k-space plane data set organized row-by-row; 

means (152) for Fourier transforming each k-space sample of the completed 
k-space plane in a first direction to recover spatial content in the first direction; and 

means (162) for Fourier transforming the k-space plane data set organized 
row-by-row in a second direction transverse to the first direction to recover spatial 
content in the second direction. 



Abstract of the disclosure 
In a magnetic resonance imaging data acquisition method, k-space (100) is 
divided into a central region (102) disposed at k-space center and one or more annular 
surrounding regions (104, 106) having increasing distances from k-space center. The 
one or more annular surrounding regions include an outermost surrounding region 

5 (106) having a largest distance from k-space center. The k-space samples in the central 
region (102) are acquired. Subsequent to the acquiring of k-space samples in the central 
region, k-space samples are acquired in the one or more annular surrounding regions 
(104, 106). The k-space samples in the outermost surrounding region (106) are 
acquired last The acquiring of k-space samples in at least the outermost surrounding 

10 region uses a row-by-row acquisition ordering in which each row of k-space samples 
acquired in the outermost surrounding region (106) completes a k-space plane. 
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